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BACKGROUND OF THE INVENTION 

Field of the Invention 

[001] The present invention relates to a multi-user detection scheme in a direct sequence 
code-division multiple-access (DS-CDMA) system and, more particularly, to a method and 
apparatus for reducing the computations of a DS-CDMA receiver. 
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Description of Related Art 

[002] With the advent of cellular telephones and portable computers, wireless 
communication has become a preferred and increasingly prevalent method of communication. 
Multi-user detection used in wireless communication is directed toward the demodulation of 
digitally modulated signals in the presence of multi-access interference (MAI). Multi-user 
detection is most often applied in Code-Division Multiple-Access (CDMA) receiver designs. 
CDMA involves a method of multiplexing a group of users where each user is identified by a 
distinct code. Multiplexing allows any number of users the ability to transmit signals 
simultaneously, and allocates the entire available frequency spectrum to each user. 

[003] DS-CDMA is a popular technique of CDMA in which the transmitter multiplies 
each user's signal by a distinct code waveform. The receiver, due to an overlap of signals in time 
and frequency, detects a signal composed of the sum of all users' signals. The MAI component 
of the received signal is the result of random time offsets between the respective users 1 signals 
and makes the optimal design of orthogonal waveforms very difficult. In detecting a signal, a 
conventional DS-CDMA system uses an auto-correlation technique in which the entire received 
signal is correlated to a particular user's code waveform. This technique is limited because of a 
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single-user detection strategy, which detects each user separately, not taking into account 
information within the signals of other users. 

[004] Multi-user detection strategies have been developed that use information of many 
| users to better detect the coded signal of each individual user. These strategies are further 

i 

described by Shimon Moshavi in "Multi-User Detection for DS-CDMA Communications", IEEE 
Communications Magazine, October 1996, hereby incorporated by reference. Some prior art 
multi-user detection methods use techniques that model the received signal over the entire 
message. While this method captures a comprehensive interaction of users, the computations 
involve an ensemble of many bits, which introduces an abundance of excess processing. Other 
prior art methods involve partial correlation, which accounts for the partial overlap of bits within 
a channel. However, such methods require much iteration to compensate for signal loss due to 
interference. 

[005] In light of the limitations of the aforementioned prior art, there is a need for a DS- 
CDMA system that can effectively reduce the effects of signal interference when detecting the 
signal of a particular user within a multi-user system. 

SUMMARY OF THE INVENTION 

[006] The present invention provides a method and apparatus to significantly reduce the 
amount of processing necessary to attain a reduction in the effects of interference when detecting 
the signal of a single-user within a multi-user system. 

[007] In accordance with the purpose of the invention, as embodied and broadly 
described, one aspect of the invention provides a method of single-user projection detection in a 
multiple access communication system. The method comprises the steps of defining a frame 
including a plurality of user signals, each user signal having at least three symbols within a 
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window of the frame; approximating a plurality of symbol amplitudes for each user among a 
plurality of users in the window; computing a set of correlation matrices based upon a plurality 
of channel parameters in the frame; calculating a current estimate of the approximated symbol 
amplitudes corresponding to at least one symbol of every user in the window using the set of 
correlation matrices in a perturbation correction technique; and shifting the window of the frame 
and returning to the step of calculating a current estimate of the symbol amplitude, when an 
entire user signal has not been processed 

[008] Another aspect of the present invention provides a method of multi-user 
decorrelation detection in a multiple-access communication system. The method comprises the 
steps of defining a frame including a plurality of user signals, each user signal having at least 
three symbols within the frame; approximating a plurality of symbol amplitudes for each user 
among a plurality of users in the frame; computing a decorrelation matrix and a set of correlation 
matrices based upon a plurality of channel parameters in the frame; calculating a current estimate 
of the approximated symbol amplitudes corresponding to at least one symbol of every user in the 
window using the decorrelation matrix and a set of correlation matrices in a perturbation 
correction calculation; and shifting the window of the frame and returning to the step of 
calculating a current estimate of the symbol amplitude, when an entire user signal has not been 
processed. 

[009] Still another aspect of the present invention provides a method of detecting at 
least one of a plurality of users in a multiple-access communication system. The method 
comprises the steps of defining a frame including a plurality of user signals, each signal being 
allocated at least three symbols within a window of the frame; comparing a delay spread of the 
frame to a predetermined threshold; approximating a symbol amplitude for each user of the 
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plurality of users in the window; computing a set of correlation matrices and a decorrelation 
matrix based upon a plurality of channel parameters in the frame, the set of correlation matrices 
and decorrelation matrix being determined by a result of the delay spread comparison; 
calculating a current estimate of the approximated symbol amplitude corresponding to at least 
one symbol of each user in the window using the determined matrices in a perturbation 
correction calculation; and shifting the window of the frame and returning to the step of 
calculating a current estimate of the symbol amplitude, when an entire user signal has not been 
i processed. 

[010] Another aspect of the present invention provides a method of detection in a 
multiple access communication system. The method comprises the steps of defining a frame 
including a plurality of user signals, each user signal having at least three symbols within a 
window of the frame; approximating a symbol amplitude for each user of the plurality of users in 
the frame; obtaining a perturbation correction value using immediate feedback, wherein a new 
estimatedampiitude is used as soon as it is available; calculating a current estimate of the 
approximated symbol amplitude corresponding to at least one symbol of every usfcr in the frame 
based upon a perturbation correction calculation; and shifting the window and returning to the 
step of calculating a current estimate of the symbol amplitude, when an entire user signal has not 
been processed. 

[Oil] In addition, another aspect of the invention provides an apparatus for detecting at 
least one of a plurality of users in a multiple-access communication system. The apparatus 
comprises means for defining a frame, wherein the frame includes a plurality of user signals each 
user signal having at least three symbols within a window of the frame; means for approximating 
a plurality of symbol amplitudes for each user among a plurality of users in the window; means 



for computing a set of correlation matrices based upon a plurality of channel parameters in the 
frame; means for calculating a current estimate of the approximated symbol amplitudes 
corresponding to at least one symbol for every user in the window using the set of correlation 
matrices and the decorrelation matrix in a perturbation correction technique; means for shifting 
the window of the frame and recalculating the current estimate of the symbol amplitude when an 
entire user signal has not been processed; and means for detecting the signal of at least one of a 
plurality of users when the symbol amplitude for the at least one user has been estimated. 

[012] Additional features and advantages of the invention will be set forth in the 
description which follows, and in part will be apparent from the description, or may be learned 
by practice of the invention. The objectives and other advantages of the invention will be 
realized and attained by the system and method particularly pointed out in the written description 
and claims hereof as well as the appended drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[013] Fig. 1 A illustrates the structure of a perturbation CDMA receiver consistent with 
the present invention; 

[014] Fig, IB illustrates the defined frame of user signals consistent with the present 
invention; 

[015] Fig. 2 shows a flow diagram illustrating single-user projection detection 
consistent with the present invention; 

[016] Fig. 3 shows a flow diagram illustrating multi-user decorrelation detection 
consistent with the present invention; 

[017] Fig. 4 shows a flow diagram illustrating a combined detection of at least one user 
consistent with the present invention; 
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[018] Fig. 5 shows a flow diagram illustrating single-user projection detection using 
immediate feedback consistent with the present invention; 

[019] Figs. 6-1 1 illustrate the results of perturbation correction simulation consistent 
with the present invention. 

DETAILED DESCRIPTION 

[020] Figure 1 A illustrates general features of a perturbation CDMA receiver 100 
consistent with the present invention. At a transmitter (not shown) a signal of a user is 
modulated using a shift-keying technique, multiplied by a distinct user code, and subsequently 
transmitted. A received signal r(t) includes information of all users transmitting information 
across multiple channels. Perturbation CDMA receiver 100 receives signal r(t) as a series of 
binary pulses called symbols. To detect a signal of one or more users, perturbation CDMA 
receiver 100 includes a user code i generator 102, user code j generator 1 16, correlator banks 
104, 1 12, 118, received signal frame buffer 108, summer banks 1 10 and 124. 

[021] Code i generator 102 is an oscillator that locally generates a distinct user code to 
detect the information of a single user / in received signal r(t) . Due to code delay derived from 
the delay lock tracking of received signal r(t) y the user code generated by user generator 102 
begins at a time T/. 

[022] Codey generator 1 1 6 is an oscillator that locally generates the distinct user codes 
corresponding to users other than user / having information contained in received signal r(t). 
User signals having distinct codes and delays are identified as interference signals. 

[023] Frame buffer 108 buffers the received signal r(t) and separates the signal into 
contiguous windows, .where each user signal within a particular window contains three symbols 
with values -1, 0, and L 
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[024] Correlator banks 104, 112, and 118 contain multiplier banks 105, 113, and 119, 
respectively for matching the information of two separate signals, and integrator banks 1 07, 1 15, 
and 121, respectively, to suppress portions of the matched information. The signals matched by 
each abovementioned correlator will be described in further detail below. 

[025] Summer banks 110 and 124 add the results of at least two signals to calculate the 
effects of interference. The operation of summer banks 1 10 and 124 will be described in further 
detail below, 

[026] In Fig. 1 A, multiplier bank 1 05 of correlator bank 104 matches a sample of the 
code generated by user code i generator 102 and a corresponding coherent sample from the 
symbols S(-l), S(0), and S(+l), respectively, of received signal r(t), through a multiply 
operation. Next, integrator bank 107 of correlator bank 104, integrates the product over the 
duration of each respective symbol, and dumps the integration result. The dumped integration 
result represents single-user detection values or zero-order amplitudes of the user signal and is 
used when calculating the amplitude of users other than user / having information contained in 
received signal r(t). 

. [027] Multiplier bank 1 19 of correlator bank 1 18 matches a sample of each generated 
user code jl through yAf generated by code j generator 1 16 corresponding to a user signal 
included in signal r(t), to a sample of the symbols S(-l), S(0), and S(+l), respectively, 
corresponding to a window of received signal r(t) > through a multiply operation. Integrator bank 
121 of correlator bank 1 1 8 integrates the product producing three zero-order amplitudes for each 
generated user code j J through jM at the given code delay. Next, multiplier bank 122 multiplies 
the zero order amplitudes by the corresponding generated user code jl through yM Next, 
summer bank 124 adds the results of multiplier bank 122 producing a signal having a time 
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sequence in parallel with the received signal rft). Summer bank 110 produces signal r 'ft) by 
computing a difference between the resultant signal of summer hank 124 and the portion of 
signal rft) contained in the window under analysis. Multiplier 113 of correlator bank 112 
matches the middle symbol of signal r'ft) with the user code i generated by user code i generator 
102 and integrator 115 suppresses any unwanted portions of the signal to produce the first order 
corrected amplitude a in . 

[028] Fig. IB illustrates an example of a 3-symbol window used in implementing the 
present invention. In practical CDMA applications, the channel is asynchronous, which means 
the non-coherent signals are randomly delayed and offset from one another. Because each user 
code has a different delay, the symbol boundary of a particular user code relative to the symbols 
of all the user codes in a window are non-uniform, i.e., asynchronous. Therefore, in 
embodiments consistent with the present invention, at time t within any 3-symbol window, the 
received non-coherent spread spectrum baseband signal corresponding to M users can be written 
as follows: 

r(D = Z^EL-.^ViCt- - nT) (1) 
where t|r,- = spreading code of user i, 
t/ = time delay of x|r, at t = 0, 
T = symbol period; t,- < T y 
din = amplitude of nth symbol of ith user, 
[029] is defined for t between x ( + n7*and i/ + (n H- 1)7J and % is zero for t outside 
the current window of the frame. As shown in Fig. IB, any three users are represented by three 
codes: Code 1, Code 2, and Code L Each user's code is segmented into three symbols (SI, S2, 
S3), and each user has one symbol starting in each of three time intervals [-T, 0], [0, T], and [T, 
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2T]. For example, user i has the first symbol SI starting at time T f - T and ending at time x f , the 
second symbol S2 starts at time t, and ends at time t,- + T, and the third symbol S3 starts at time 
it + T and ends at time T, + 2T. Any symbol that does not fall in current intervals will have no 
effect on the symbols in the current intervals. 

[030] A detailed flowchart illustrating one aspect of perturbation correction calculations 
according to the present invention is shown in Fig. 2. The flow chart illustrates the processing 
that takes place during single-user projection by the CDMA receiver 100. With reference to Fig. 
2, in step 200, frame buffer 108 receives die transmitted signal and defines the window as shown 
in Fig, IB. Over the course of the transmission, the user signals within a current window of the 
frame are communicated over a significantly invariant channel. Specifically, the channel 
parameters for each user including code, code delay, carrier frequency, and phase are nearly 
constant over the entire signal. Code delay is the beginning time of a code measured by a 
common reference clock. Carrier frequency is the frequency of the sinusoidal function of the 
oscillator (not shown) used in the transmitter to up convert the baseband signal to a desired band 
(range of frequency) for transmission. CDMA receiver 100 regains the original baseband signal 
by down converting the received signal with a signal generated by user code / generator 102 
having the same carrier frequency and phase of the original oscillator. The variable r in Eq. (1) 
defines the baseband signal. The correlation value is an implicit function of the measured r(t) 
and is usually a complex quantity due to the variable phase at the CDMA receiver 100. Here, 
although computations using r(t) mainly involve the real part of the value, the imaginary part 
may also be treated in the same manner. Therefore, embodiments consistent with the present 
invention are constructed based on the assumption that the channel parameters are constant over 
an entire frame. Each channel parameter is further defined by a set of margins — a range of 
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values having upper and lower Hmits-that also determine the duration of the frame. Because the 
values in the range are usually determined by the resolution of the measurement apparatus, the 
range values are regarded as nearly constant. The tracking of these parameters has been treated 
extensively in CDMA systems, is well known, and is not discussed in further detail. 

[031] In step 202, the amplitude of one of the user signals within the defined frame is 
approximated using an auto-correlation technique. This technique is effective to maximize the 
signal and minimize the interference or cross-correlation. The zero-order approximation is 
performed neglecting interference, which means that delay tracking is assumed perfect and user 
codes jl through 7 Mare nearly orthogonal. The single user detector computes symbol amplitude 
din with bit (± 1) data from equation (1) as follows; 



( Vi (r i + nT), Vi (T i + nT)) ;ri = C/0 
[032] Eq, (2) uses the scalar product notation (x,y) for correlation defined as (x,y) = 

eT 

x(t)y{t)dt. These terms in equation (2) are defined above in connection with equation (1). 

[033] Because this approximation is calculated assuming orthogonal codes, which are 
not normally encountered in actual practice, a correction value is obtained according to the 
perturbation theory. Using perturbation theory and assigning c'^ as the correction value and c in 
as the estimated value given by Eq. (2), a in is defined as follows: 



(2) 



Q 'm Cin ^~ C in 



(3) 



[034] If ai n defined in equation (3) is substituted into equation (1) the following equation 
(4) results: 



IlXl-^M- t,- nT)=rf0-E" E^c^Ct- t r nT) 



(4) 
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[035] In equation (4), r(t) represents the measured value of the received signal The 
variable c& represents the estimated value obtained from Eq. (2), and ty\ represents the given user 
codes. All values on the right hand side of equation (4) are known, whereas, c'^ values located 
on the left hand side of equation (4) are the only unknown values. 

[036] The perturbation correction value c'*« is calculated (step 204) based upon the 
cross-correlation between the symbol of one user and the adjacent symbols of other users within 
a current window of the frame. This correction value is obtained by calculating c',> in equation 
(4) using the present values for user /, and current symbol n = 0. When the value of the right- 
hand side of equation (4) has been calculated for the symbol of a user 7, a value of the 
perturbation correction in matrix-vector notation is approximately given as: 



c 'jo 88 tt— IwEU^HffCiOjy = ; to 



M 



(5) 



where the correlation matrix Kfi(n) represents the cross-correlation between the adjacent symbols 
n = -1, 0, 1 within the current window of the frame and between user code i and user code y. The 
correlation matrix is given by 



(6) 



where ty 0 provides a normalization value defined by N jt> = (r y )) . As 

indicated in the formula, the correlation matrix is a function of codes i and j and their respective 
delays x t and Xj and it is independent of other channel parameters. 

[037] In Fig. IB for example, the middle symbol of user 2 overlaps with other users 
with respect to at most two symbols in the frame. When the two adjacent user symbols are 
known, the right-hand side of equation (4) allows an explicit calculation of cross-correlation to 
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the middle symbol of user 2. The right-hand side of equation (4) includes r(t), the total received 
signal of all users and c,„, the estimate on the first approximation of the user's symbol amplitude. 

[038] Once the correction value has been determined by the above projection method, 
processing continues at step 206 where an estimate of the symbol amplitude is calculated. This 
calculation is performed by substituting the correction value obtained at step 204 using equation 
(5) for the term c' in in equation (3). At step 208, this estimate of symbol amplitude replaces 
the first approximation of the symbol amplitude obtained in step 202. 

[039] At step 210, it is determined if the entire signal of the current window has been 
processed, and every user's amplitude has been corrected. When the entire signal has been 
processed, i.e., the amplitude of the entire signal has been estimated, the final estimated 
amplitude value is obtained from equation (2), the value of a m allows the approximate value 
of r(t) to be obtained. In contrast, if the entire signal has not been processed, processing 
continues at step 212 where it is determined whether a new frame has been defined (or obtained, 
or present, or received). When a new frame is defined, processing returns to step 202 where a 
first approximation of the signal amplitude is performed and steps 204 through 210 are repeated. 
However, when a new frame is not defined (step 212), the next iteration of amplitude estimation 
is performed. This next iteration updates the estimated amplitude value based upon a new 
symbol amplitude estimate of the current user in the current window of the frame. Performing 
the calculation of the next iteration begins at step 214 by shifting the current window by one user 
symbol. This window shift is defined in equation (5) by the term n that corresponds to samples 
of user code i. Processing continues by returning to step 204, where the correction value 
corresponding to the new user symbol within the current window of the frame is calculated and 
steps 206 through 210 are repeated. 
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[040] A second aspect of the present invention is illustrated by the flow chart in Fig. 3. 
The flow chart of Fig. 3 illustrates the processing of the CDMA receiver 100 during multi-user 
detection. The processing in steps 300, 302, and 306 through 314 are consistent with the 
processing as described above in steps 200, 202, and 206 through 214 of Fig. 2, respectively. At 
step 304, a decorrelation matrix is used to obtain the perturbation correction value. Perturbative 
decorrelation is performed by including the correlation of terms Y/Ty) with the middle term 
(«=0) on the left-hand side of equation (4) to yield the following equation: 

where c) Q is obtained by inverting the matrix HrfQ). The inversion matrix (Hji(0)f f is 
representative of the decorrelation matrix and Hj,(n) is representative of the correlation matrices. 
Use of the decorrelation matrix is most effective when the signals of the users are nearly 
synchronous. As noted above, there is no overlap of user signals within a synchronous system, 
i.e., the signals of other users contribute no useful information to each other. Therefore, for 
simplicity and limitation of M equations for M unknowns, the terms n-A and n- 1 in the sum 
of the left-hand side of equation (4) are neglected. Moreover, when n=0 the correlation matrix 
Hji(0) is a subset of correlation matrix Hji(n). 

[041] The flow chart of Fig. 4, illustrates another process consistent with the present 
invention for performing high order perturbation calculations. The initial steps (400, 402) of 
defining a frame and approximating a symbol amplitude for a particular user are performed as 
explained above regarding the flowcharts in Figs. 2 and 3. Next, at step 404, a delay spread of 
the symbols within the current window of the frame is evaluated to determine the type of matrix 
to be used during perturbation correction. If the delay spread is small, e.g., less than two chips, 
there is very little overlap between dropped and included terms, therefore multi-user 
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decorrelation detection is performed according to steps 418-428, which correspond to previously 
described steps 304-314, respectively, using a correlation matrix. Otherwise, single-user 
projection detection according to steps 406-416, which correspond to previously described steps 
204-214, respectively, is performed using a decorrelation matrix. The multi-user decorrelation 
and projection detection schemes operate as previously described for the process shown in Figs. 
2 and 3, However, when a new frame is defined, operation in each processing scheme returns to 
step 402 where an approximation of the symbol for the current user is acquired. Step 404 is 
again performed when a new frame is introduced. 

[042] The flow chart shown in Fig. 5 illustrates a further process consistent with the 
invention. The process in Fig. 5 relates to the perturbation method of Fig. 2, in the case in which 

u 

j| the value of the previous approximation is used to calculate the next symbol. Steps 500, 502, 
and 506-514 are the same as steps 200, 202, and 206-214, respectively, of Fig. 2 and 300, 302, 
and 306-314 respectively of Fig. 3. In Fig. 5 at step 504, a perturbation correction value is 
obtained using an immediate feedback technique. This feedback technique is defined by 
equation (8) where the summation on the right-hand side (RHS) of equations (5) and (6) is 
changed to the following expression: 

Z( )- ImILo^^+Im^..) 11 /- 1 ) («) 

Using immediate feedback for the calculation of the next symbol significantly improves the 
perturbation correction on the basis that there is no delay about the correction. This is in contrast 
to the prior art where many symbols are processed before feedback can be performed* 

[043] For the processing schemes described above and illustrated in Figs. 2-5, obtaining 
the correlation or decorrelation matrix for a whole frame requires substantial computation. The 
computations involve all rows and columns of all users for each matrix. As noted above, when 
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the sigaal codes of users are orthogonal and are identical for each symbol, the number of 
computations is relatively small. However, if each symbol of the orthogonal code is multiplied 
by a scrambling code, which is longer than the symbol and overlaps adjacent frames, then a 
matrix of every symbol must be computed. A scrambling code is a random sequence multiplied 
chip by chick on the spreading code. This is done to further randomize the spreading codes so 
that they have less cross correlation when they are not synchronous. In the multi-user detection 
receiver, scrambling becomes unnecessary because the interference will eventually be cancelled. 
This procedure invariably produces a very large number of computations. On the other hand, 
when the scrambling code is longer than the symbol but much shorter than the frame, the 
calculation of the correlation matrices is still efficient because the same correlation matrix can be 
reused. 

[044] In accordance with the present invention, in order to address the computational 
complexity resulting from use of a longer scrambling code, instead of calculating the correlation 
matrix for every symbol, it is more efficient to solve equation (4) directly because no correlation 
matrix is needed. The summation on the right-hand side of the equation is performed sample by 
sample, as shown in Fig. 1 A, and both sides are correlated with fyfa ) to obtain cp. 

[045] Fig. 6 illustrates the simulation of a static model in which each user transmits a 
code at a fixed and equal absolute amplitude and randomizes his delay at the beginning of a 
frame. The static model shows 16-dimensional Walsh codes being used by 1 6 users, and each 
user communicating 100 bits over the system. The 100 bits are divided into 10 frames, each : 
consisting of 10 bits. Therefore, the signal of each user is measured in 10 bits/frame increments 
until the signal over the entire frame has been simulated. 
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[046] Figs. 7 and 8 illustrate the results of simulating a near/far model using the present 
invention. Near/far is simulated by identical independent distribution (iid) of the user's 
amplitude, a technique well known in the art. All channel parameters including code delay and 
absolute amplitude of a user stay constant in a 10-bit frame but changes from frame-to-frame. 
The model illustrates the near-far problem, where signals of users transmitted from a location 
close to a receiver undergo less attenuation than those farther away. As a result, the signals 
transmitted from a farther distance encounter greater interference, especially from the signals 
transmitted by closer users. Figs. 7 and 8 illustrate near-far resistance for the case of 16 and 32 
users, respectively. 

[047] Higher order perturbation processing can further improve results. In high order 
perturbation processing, the current estimated value c« in Eq. (4) is equal to the symbol 
amplitude approximation value <2& of the lower order perturbation. As discussed previously, Eq. 
(2) provides the amplitudes resulting from zero order pertuibation. Eqs. (5) and (7) provide the 
first order perturbation based on projection and decorrelation and also provide identical formulas 
for higher order perturbation. In the second order perturbation, the estimated value c in> in the 
summation of Eqs. (5) and (7), represents the amplitude of the first order perturbation. 
Regarding the first order perturbation, the first term rj in Eqs. (5) and (7) is the same term 
correlated with the middle symbol in the 3-symbol window during single user detection. While 
Fig. 1 A shows only the first order perturbation, Figs. 2-5 illustrate any order of perturbation in 
terms of the number of iterations executed in a frame before a new frame is triggered. If the 
results of the perturbation converge, i.e., c \ n approaches zero, then the number of iterations is 
limited. Particularly, the convergence is guaranteed when the signal to interference ratio is larger 
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than a predetermined value that is greater than one. High order perturbation processing 
determines the next (higher) order of symbol amplitude approximation, using Eqs. (3)-(8). 

[048] Fig. 9 is a snapshot of Fig. 8 at a point in time in which the delay spread of the 
signals is small. 

[049] Figs. 10 and 1 1 show the results of higher order perturbation simulation for 16 
users. Particularly, Fig. 10 shows the results using a projection method, whereas Fig. 1 1 shows 
the results using decorrelation. Additionally, the third order perturbation results illustrated in 
each figure exhibit significant improvements over the first order perturbation results. 

[050] Equations (1) through (8) as mentioned above can be implemented in hardware 
(e.g., ASIC) or software (e.g., DSP) whichever is more convenient. 

[051] It will be apparent to those skilled in the art that various modifications and 
variations can be made in the apparatus and method of the present invention without departing 
from the spirit or scope of the invention. Thus, it is intended that the present invention cover the 
modifications and variations of this invention provided they come within the scope of the 
appended claims and their equivalents. 
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